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Abstract 

In this paper, the author present reliable symbolic algorithms for solving a general bordered tridi- 
agonal linear system. The first algorithm is based on the LU decomposition of the coefficient matrix 
and the computational cost of it is 0(n). The second is based on The Sherman-Morrison- Woodbury 
formula. The algorithms are implementable to the Computer Algebra System (CAS) such as MAPLE, 
MATLAB and MATHEMATICA. Three examples are presented for the sake of illustration. 
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1 Introduction 

A general bordered tridiagonal linear system takes the form: 
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(1.1) 



Many problems in mathematics and applied science require the solution of a general bordered tridiagonal 
linear system. For example the solution of certain partial differential equations, spline approximation, 
computation of electric power system, etc. [1-4]. Recently in [5], the author presents an approach to find 
solution of the equation (1.1) in heterogeneous environments. In this approach, a bordered tridiagonal 
linear system is first converted into three or more tridiagonal linear systems that are independent each 
other, then the solution of the bordered tridiagonal linear system is obtained by solving the tridiagonal 
linear systems via parallel computing under heterogeneous environments. The motivation of the current 
paper is to establish efficient algorithms for solving a general bordered tridiagonal linear system of the form 
(1.1). 

The paper is organized as follows: In Section 2, a new symbolic algorithm is constructed. The Sherman- 
Morrison- Woodbury algorithm is presented in Section 3. Three illustrative examples are given in Section 
4. Conclusions of the work are presented in Section 5. 
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2 A new Symbolic Algorithm 



In this section we shall focus on the construction of new symbolic algorithm for solving a general 
bordered tridiagonal linear system. To do this we begin by considering the LU decomposition [6] of the 
coefficient matrix in (1.1) 
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where di, and /3i in (2.1) satisfy: 
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We also have the determinant of the coefficient matrix in (1.1): 

n 

Determinant = di . 
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At this point it is convenient to formulate our first result. It is a symbolic algorithm for computing the 
solution of a general bordered tridiagonal linear system (1.1). 

Algorithm 2.1 the solution of a general bordered tridiagonal linear system (1.1), we may proceed as fol- 
lows: 

INPUT number of equations in (1.1) n and the components ai, i = 1,2, n, bi, i = 1, 2, n — 1, 

Ci, i = 2, 3, n, pi and q i} i = 1,2, n — 2. 
OUTPUT The solution of a general bordered tridiagonal linear system x. 

Step 1: Set di = ai, ol\ = j-, Pi = pi. If di = then d\ =t(t is just a symbolic name) end if. 
Step 2: Fori = 2,3,..., n- 1 

Compute di = ai — * Cj/rfj_i ; 
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If di = then di = t end if. 
Step 3: For i = 2,3, n — 2, compute 

a% = (qi - oti-i * bi-i)/di, 
Pi=Pi- Pi-i * Ci/di-i. 
Step 4: Set a n -i = (c„ - a„_ 2 * b n - 2 )/d n -i, 
Set B n -i = 6„_i - /3„_ 2 * c n -i/d n - 2 , 
Set d n = a n — Y^j=i a j * fij ■ If d n = then d n — t end if. 
Step 5: Compute Determinant = ^n™=i dij 
Step 6: Set z\ = yi . 
Step 7: For i = 2, 3, n — 1, compute 

Zi=D%- Ci* Zi-i/di-i- 
Step 8: Compute z n = y n - YTjZi a i * z n x n = z n /d n , and x n -\ = (z n -i - f3 n - 1 x n )/d n - 1 . 
Step 9: For i = n — 2, n — 3, 1, compute 

Xi = (zi - b t x l+1 - j3 i x n )/d l . 
Step 10: Substitute t = in all expressions of the solution vector x. 

The symbolic Algorithm 2.1 will be referred to as SBTLS algorithm. The computational cost of 
SBTLS algorithm is 19n — 34 operations. In [7], the PERTRI algorithm is special case of our algorithm 
when pi = qi = 0, i = 2, 3, • • ■ , n — 2. 



3 The Sherman-Morrison- Woodbury algorithm 

In this section, we are going to formulate a new symbolic algorithm for solving a general bordered 
tridiagonal linear system of the form (1.1) based on the Sherman-Morrison- Woodbury formula and any 
symbolic tridigonal linear solver. 

A general bordered tridiagonal linear system of the form (1.1) can be written in the form: 
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Thus (3.1) is equivalent to: 



Mix +Vx = y 

„ „ (3-2) 
U 1 x + M 2 x = y 

Assume that a n ^ 0. After elimination of x from (3.2), we get the linear system: 

Mx = y (3.3) 



■A 



where 

M = Mi- V M 2 " 1 U T , and y = y - VM 2 ~ 1 y" . 

If we applying the Sherman-Morrison- Woodbury formula[8] to M, we will obtain: 

M- 1 = Mf 1 + Mf V (M 2 - £/ T Mf V)" 1 U T M^\ 
and 

a;' = M~ x y = r + Mf V (M 2 - J7 T Mf V)" 1 C/ T r. 

where r is the solution of Mir = y. It is clear that the solution x can be found from the above formula 
by successive calculation of the expressions 

r = MfV, 9 = Mf 1 V r , £/ T MfV, (M 2 - f7 T Mf V)" 1 , and (M 2 - t/ T Mf V)" 1 C/ T r. 

The main part of above calculations is finding the first two expressions, which is equivalent to solving 
two (n - l)-by-(n - 1) tridiagonal linear systems with the same coefficient matrix Mi and different right- 
hand sides. After finding of x , we can get x from the second equation of (3.2) by formula 

x" = Mf 1 {y" - U T x'^j 



At this point it is convenient to formulate our second result. It is a symbolic algorithm for computing 
the solution of a general bordered tridiagonal linear system of the form (1.1) and can be considered it as 
a natural generalization of the The Sherman-Morrison algorithm in [7] . 



Algorithm 3.1 To compute the solution of a general bordered tridiagonal linear system (1.1), we may 

proceed as follows: 
Step 1 : Find M 1 ,M 2 ,U T ,V,y ,y" , and y = y - VM% 1 y " . 
Step 2: Solve M x r = y, and M x q = V. 

Step 3: Compute x = r + q (M 2 - U T q)^ 1 U T r, and x" = M^ 1 (y" - U T x'^j . 

X 

Step 4 : Compute the solution x — ,, 



Two systems Mir = y, and M\q = V in algorithm 3.1 can be solved in parallel by any symbolic 
tridiagonal algorithm. The symbolic Algorithm 3.1 will be referred to as SMWBTLS algorithm. 



4 ILLUSTRATIVE EXAMPLES 

In this section we give three examples for the sake of illustration. 
Example 3.1. Let 
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■1.0861, 2.6376, -3.0065] 5 



i) By applying the BTLE algorithm in [5] , it yields 

xbtle = [3.8709,-2.2703,3.1714,1.8793, 

ii) By applying the SBTLE algorithm, it yields 

xsbtle = [3.8638, -2.2838, 3.1464, 1.9121, -1.0871, 2.6192, -2.9767] T . 

iii) By applying MATLAB command A\b, it yields 



A 



xmatlab = [3.8638, -2.2838, 3.1464, 1.9121, -1.0871, 2.6192, -2.9767] T . 
iv) By applying SMWBTLE algorithm , it yields 

xsmw = [3.8638, -2.2838, 3.1464, 1.9121, -1.0871, 2.6192, -2.9767] T . 



Example 3.2. Let 
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i) By applying the BTLE algorithm in [5] , it breaks down 

ii) By applying the SBTLE algorithm , it yields 
xsbtle 
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[1,1,1,1,1,1,1,1,1,1] T . 



Example 3.3. We consider the following general bordered linear system in order to demonstrate the 
efficiency of SBTLE algorithm and SMWBTLE algorithm. Let 
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It can be verified that the exact solution is x = [1, 1, • • • , 1] T . We used proposed algorimths, SBTLE and 
SMWBTLE Algorithms, and Gauss Elimination Algorithm to compute the solution, x. Results are given 
in the next table in which e = llx — xlloo- 



Table!. 



n 


e = \\x — x |oo and CPU time 


SBTLE Algorithm 


SMWBTLE Algorithm 


Gauss Elimination Algorithm 


500 


3.4100 x 10-«j 0.0613 


1.4600 x 10 7 | 0.1633 


4.4000 x 10~ 8 | 0.0616 


000 


6.9100 x 10" 8 0.1113 


2.9600 xl0~ Y 0.3133 


8.9000 x 10~ 8 0.1116 


5000 
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4.4900 x 10~ 7 0.5116 


10000 


6.9910 xl0- y | 1.0113 


2.9960 x 10- b 3.0133 


8.9900 x 10-' 1-0116 



5 CONCLUSIONS 

In this work new symbolic computational algorithms have been developed for computing the solution of 
a general bordered tridiagonal system. The algorithms are reliable, computationally efficient and remove 
the cases where the numeric algorithms are fail. 
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